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ABSTRACT 

Radial-velocity measurements and sine-curve fits to the orbital radial velocity vari- 
ations are presented for ten close binary systems: TZ Boo, VW Boo, EL Boo, VZ CVn, 
GK Cep, RW Com, V2610 Oph, V1387 Ori, AU Ser, and FT UMa. Our spectroscopy 
revealed two quadruple systems, TZ Boo and V2610 Oph, while three stars showing 
small photometric amplitudes, EL Boo, V1387 Ori, and FT UMa, were found to be 
triple systems. GK Cep is close binary with a faint third component. 

While most of the studied eclipsing systems are contact binaries, VZ CVn and 
GK Cep are detached or semi-detached double-lined binaries, and EL Boo, V1387 Ori 
and FT UMa are close binaries of uncertain binary type. The large fraction of triple 
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and quadruple systems found in this sample supports the hypothesis of formation of 
close binaries in multiple stellar systems; it also demonstrates that low photometric 
amplitude binaries are a fertile ground for further discoveries of multiple systems. 

Subject headings: stars: close binaries - stars: eclipsing binaries - stars: variable stars 



1. INTRODUCTION 

This paper is a continuation of a series of papers (Papers I - XIII) of radial- velocity studies of 
close binary stars and presents data for the thirteenth group of ten close binary stars observed at 
the David Dunlap Observatory (DDO). Because of the closure of the observatory, it is most likely 
the last paper which retains the usual format of 10 new orbits per paper; the last paper, Paper 
XV, will conclude the series with results for partly covered systems and for variable stars which 
were found no t to b e binaries. For full references to the previous papers, see the last paper by 



Rucinski et al.l (|2008l . Paper XIII); for technical details and conventions, for preliminary estimates 



of uncertainties, and for a descr iption of the broadening functions (BFs) technique, see the interim 



sumrnary p aper iRucinskil (|2002l . Paper VII). The DDO studies have used the efficient program of 



PvchI (|2004l ) for removal of cosmic rays from 2-D images. 



All data used in the present paper were obtained using the broadening functions extracted 
from the region of the Mg I triplet at 5184 A, as in most of the previous papers, using the new 2160 
lines/mm grating acquired at DDO in August 2005. The radial- velocity (hereafter RV) observations 
reported in this paper have been collected between May 2006 and the memorable day of July 2, 
2008 when the David Dunlap Observatory ceased to operate. The ranges of dates for individual 
systems can be found in Table [U for TZ Boo, where we used the smoothed BFs, the dates are in 
Table [3] giving RVs for the companion of the short-period binary. 

Throughout our program, selection of the targets was quasi-random: At a given time, we 
observed a few dozen close binary systems with periods usually shorter than one day, brighter 
than 10 - 11 magnitude and with declinations > —20°; we published the results in groups of ten 
systems as soon as reasonable orbital elements were obtained from measurements evenly distributed 
in orbital phases. 

Among the present targets, three systems, AU Ser, VW Boo and VZ CVn, have had reliable 
radial velocity orbits previously published. In addition, EL Boo and FT UMa were originally 
classified as pulsating variables while V1387 Ori does not have any ground-based photometry, but 
only Hipparcos data. More details are given in Section [2] in the descriptions of the individual stars. 

The RVs for the short period binaries reported in this paper were determined by fitting the 
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double rotational profiles to the extracted BFs, as explained in iPribulla et al.l (|2006l ). Similarly 
as in our previous papers dealing with multiple systems (here the six cases of TZ Boo, EL Boo, 
GK Cep, V2610 Oph, V1387 Ori, and FT UMa), the RVs for the eclipsing p air were obtained afte r 
removal of the slowly rotating components, as was described most recently in lRucinski et al.l (|2008l ). 
In the case of GK Cep, where the third component is faint and not well visible in all spectra, its 
average profile was subtracted before the binary star analysis. 

As in other papers of this series, whenever possible, we estimated spectral types of the program 
stars using new classification spectra centered at 4200 or 4400 A. Good quality, multiple classifica- 
tion spectra were obtained before the closure of the observatory only for EL Boo, VZ CVn, GK Cep, 
V1387 Ori, and FT UMa; for the remaining targets the classification spectra were obtained only 
once and in poorer conditions so that the spectral types are less reliable. The spectral types were 
comp ared with the mean (B — V) color indices usually taken from the Tycho-2 catalog (IHeig et al 



2000 ) and the photometric estimates of the spectral types using the relations of Besselj ( 19791 ). In 
this paper we also rn ade use of infrared colors determined from the 2^ All Sky Survey (2MASS, 
Skrutskie et al.l 120061 ). Especially useful is th e J — K co lor index, which is monotonically rising 



from the early spectral types to about MOV (|Goxl |2000| ). This infrared color is less affected by 
the interstellar absorption than the B — V index. Paralla xes cited througho ut the paper were 
adopted from the new r eduction o f the Hipparcos raw data (jvan LeeuwenI 120071 ) which supersedes 
the original reductions (lESAlll997l ). 



This paper is structured in a way similar to that of previous papers, in that most of the data 
for the observed binaries are in two tables consisting of the RV measurements in Table [Hand of their 
preliminary sine-curve solutions in Table [2j Radial velocities and the corresponding spectroscopic 
orbits for all ten systems are shown in phase diagrams in Figures [H - El The measured RVs are 
listed in Table [TJ Table [2] contains also our new spectral classifications of the program objects. 
Section [2] of the paper contains summaries of previous studies for individual systems and comments 
on the new data. Examples of BFs of individual systems extracted from spectra observed close to 
quadratures are shown in Fig. [H 

The data in Table [2] are organized in the same manner as in the previous papers of this series. 
In addition to the parameters of spectroscopic orbits, the table provides information about the 
relation between the spectroscopically observed upper conjunction of the more massive component. 
To (not necessarily the primary eclipse) and the recent photometric determinations of the primary 
minimum in the form of the O — C deviations for the number of elapsed periods E. The reference 
ephemerides were taken from various sources: for EL Boo, and FT UMa, we doubled the Hipparcos 
period and shifted the instant of the maximum by 0.25P : for V1387 Ori we us ed the Hipparcos 
ephemeris; and for V2610 Oph, it has been adopted from lWils &: Dworald (120031 ). For the rest of 
the systems, th e ephemerides given in the on-line version of "An Atlas of O-G diagrams of eclipsing 
binary stars"EI (Kreiner 2004) were adopted. As the on-line ephemerides are frequently updated. 
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we give those used for the computation of the O — C residuals below Table [2] (up-to-date as of July 
2008). The deeper eclipse in W-type contact binary systems corresponds to the lower conjunction 
of the more massive component; in such cases the epoch in Table [2] is a half-integer number. 



2. RESULTS FOR INDIVIDUAL SYSTEMS 



2.1. TZ Boo 



TZ Boo (HIP 74061, BD-|-40°2857) is a well known contact binary discovered bv lGuthnick &: Prager 



(|l926l ). The system is unusual: Its depths of the minima change, switching betw 'een the A and W - 
type light curve (LC) types and the LC shows very large variations of its shape (|Hoffmannlll978al ). 
In spite of total eclipses, the large variability of the LC has made a reliable photometric solution 
of the system impossible; thus, there exists no modern LC synthesis solution assuming the Roche 
model. The orbital period of the system is also variable. The cont inuous orbital period decre ase 
of TZ Boo, accompanied by wave-like variations was interpreted by iPribulla &: Rucinskil (j2006l ) in 
terms of light-time effect caused by a third component on a 34 years orbit. But an adaptive optics 
search for close visual companions to contact binaries (jRucinski et al.l 120071 ) did not show any close 
companion to TZ Boo. 



TZ Boo is a very difficult system for spectroscopic studies. IChangl (| 19481 ) found that spectral 
lines do not double in quadratures, but noticed that the measured RV of the system varies in 
the range of about 60 km s~^. The author did not fiii d any c orrelation between the RV and 
the predicted photometric orbital phase. Later iHofFmannI (jl978bl ) tried to interpret th e results of 
Chang (j 19481 ) by systemic velocity changes induced by a third body on a 10.31 day orbit. iHoffmann 
(jl978bl ) added a new set of spectroscopic observations to the older observations, but could not solve 
the problem. The author estimated the mass ratio as q = 0.2 (which in fact is identical to our 
current result) and inte rpreted the inexplic a ble b ehavior of the profiles by the enhanced activity 
in the system. Finally, iMcLean &: HilditchI (119831 ) obtained the first, seemingly fully satisfactory 
spectroscopic orbit of the system: Vq = —36.7 km s^^, i^i = 33 ± 7 km s^^, and K2 = 249 ± 38 
km s~^ resulting in g = 0.13 it 0.03. The authors did not notice a third component in the profile 
that - as we suspect now - was blended with the primary component; unfortunately this means 
that this spectroscopic orbital solution is of limited use. 

Our spectroscopy shows that TZ Boo is either a triple system with a third body on a 9.48 day 
orbit or - more likely - a quadruple system consisting of a contact eclipsing binary and a second, 
non-eclipsing, single-line binary with the period of 9.48 days. The system was found to be fairly 
difficult for the DDO telescope due to the short orbital period of the contact binary, P12 = 0.297 
days, its relatively low brightness, Vmax = 10.4, and the obviously composite spectra. To determine 
reliable parameters, we took as many as 215 spectra spanning more than two years. First, we fitted 
a three-Gaussian model to each BF and subtracted the third component contribution. The BFs of 
the contact binary were noisy so we decided to smooth and rebin them with a 0.01 step in orbital 
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phase. The resulting mass ratio, q = 0.207 it 0.005, is rather smah and inconsistent with the large 
photometric amplitude of the system (note however the disparity of the values: 0.59 mag in the 
General Catalogue of Variable Stars, 0.39 mag in Hipparcos photometry, and 0.35 mag according to 



Fig. 1 of lSchaublll990l ). It is possible that the pronounced LC changes and the W/A-type switching 
are caused by the activity and/or by relatively slow (and of unexplained nature) pulsations of the 
third component. Also, the possibility of eclipses in the second bina r y can not be ruled out. The 



changes of the LC amplitude - as seen best in Fig. 2 of lAwadalla et al.l (|2006l ) - could also be caused 
by significant long-term brightness changes of the third/fourth component. Another explanation 
could be the precession of the eclipsing-pair orbit, but this is less likely because all the observed 
LCs have shown total eclipses. The system deserves a dedicated study of the light curve changes in 
a standardized photometric system to relate the relative LC shape changes to its overall brightness. 

The orbit of the second (or outer) binary is circular - which is rather surprising considering 
the long orbital period - and is very well defined with = 43.15 it 0.14 km s~^ (see Fig. [5] and 
Tabled!). The center-of-mass velocities of the contact binary, V^'^ = —46.57 it 0.90 km s~^, and the 
non-eclipsing binary, Vq^ = —54.64 it 0.12 km s~^ are different, indicating a slow orbital motion. 
The light contribution of the third component around the quadratures of the contact binary is 
L^/{Li + L2) = 0.28 ± 0.05. Because of the well recognized problems with a proper continuum 



rectif ication of the spectra of combined broad and sharp line components (jRucinski &: Pribulla 



20081 ). we regard this estimate of the third light as an upper limit. 

The RVs of the contact pair show a rather large point-to-point scatter of about 10 km s~^. The 
scatter does not decrease under the assumption that TZ Boo is a tight triple system and that the 
contact binary revolves around the common center of mass with the third component. Although 
with the present data we cannot fully exclude this possibility, it is much more probable that TZ Boo 
is a hierarchic quadruple consisting of two binaries revolving in 34-year orbit, as indicated by the 
cyclic period changes. 

2MASS infrared color, J — K = 0.454 corresponds to G7V spectral type. Our single spectral 
classification spectrum indicates a range of admissible spectral types between late F and G5. The 
Hipparcos parallax, tt = 6.63 it 1.54 mas is unfortunately too uncertain to be of value. 



2.2. VW Boo 



VW Boo (HI P 69826, CSC 908 1 170) is a short-pe riod (P = 0.3422 days) close binary that 
was discovered bv iHoffmeisteii (j 19351 ). iBinnendijkl (jl973l ) obtained and analyzed BV photoelectric 
photometry of the system, which shows that components have rather different temperatures for 
such a close binary. Therefore, VW Boo belongs to a small group of short-period contact binaries 
that have poor thermal contact (others are e.g., CN And, FT Lup, V432 Per or AU Ser). 



Rainger et al.l (|l99d ) re-analyzed LCs of lBinnendiikI (jl973l ) and obtained the first spectroscopic 



data for the system. The authors interpreted the LCs by a contact configuration with a hot spot 
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on the secondary component, about 640 K hotter than the surrounding photosphere. The CCF 
analysis of the spectra yielded the following spectroscopic elements for VW Boo: Ki = 99.2 it 2.1 
km s-\ K2 = 230.1 ± 5.4 km s'^ (resulting in g = 0.428 ± 0.030), V^i = 21.5 ± 1.5 km s'^ and 
V m = 26.3 ±4.2 k m (determined separately for the two components). Independent observations 



of lHrivnaM (jl993l ) resulted in q = 0.45 and the total projected mass of (Ml + M2) sin^ i = 1.34 Mq. 



Our spectroscopic elements (Table [2]) are in a good accord with the results of iRainger et al 



(jl990l ). However, the extracted BFs do not show any irregularities which could indicate the presence 
of either hot or cool spots on the surface. The rotational profiles of the components in BFs taken at 
the orbital quadratures suggest that system is either in marginal physical contact or is a detached 



one. 



Although the system was observed by Hipparcos, its trigonometric parallax, vr = 2.12 it 2.52 
mas is too uncertain to be of any use. The 2MASS infrared color of the system, J — K = 0.467 
corresponds to the G8V spectral type which agrees with our classification G5V. 



2.3. EL Boo 



Variability of EL Boo (HIP 72391, SAO 101223) was detected during the Hipparcos mission, 
where it was classified as a (5 Set variable with P = 0.2068860 days. Later, it was identified as 



literature (Senavci et al 



2007 



an X-ray source (|Mickaelian et al.l 120061). Observations of only two minima are available in the 



Hubscher k Walter! [2OO7I ). In the AS AS databasfll it is classified 



as 

an eclipsing binary. The ASAS light-curve indicates a W UMa-type classification and total eclipses. 
No other photometry or spectroscopy of the target has been published yet. 

The BFs (see Fig. clearly show that EL Boo is a triple system containing a very close 
(possibly contact) binary and a slowly-rotating third star. The system was not previously known 
to be a visual binary. The orbital period of the eclipsing pair, P = 0.413772, is two times longer 
than the previous Hipparcos (pulsation) period. The third component contributes about a half 
of hght to the total brightness of the system, ^3/(^1 + L2) = 1.00 ± 0.08. The RV of the third 
component was found to be constant, RV3 = —22.6 ± 1.1 km s^^, and rather close to the systemic 
velocity of the close eclipsing binary, Vq = —24.6 it 1.1 km s~^ indicating a physical bond. The 
projected rotation velocity of the third component is low but detectable even at our rather moderate 
spectral resolution, vs'mi ~ 23 km s^^. The outer orbital period of this triple must be at least 
several years or longer. 

The relatively small photometric amplitude of the system, 0.19 mag, results from the significant 
light contribution of the third component and from the rather low mass ratio of the binary, q = 
0.248 ± 0.007. On the other hand, the projected mass of the close binary, (Mi + M2)sin^z = 
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1.448 lb 0.028 Mq is high for the F5V spectral type indicating a high inchnation angle and the 
implied possibility of total eclipses. The photometric amplitude of the close binary, without the 
diluting effect of the third light, is expected to be about 0.42 mag, which is unexpectedly large, 
even if eclipses of the close binary were total; this discrepancy, which is similar to that in TZ Boo, 
indicates that we may have overestimated the amount of the third light. The system definitely 
deserves new, high-precision photometric observations. 

The published F8 spectral type is fairly late for a contact binary with a 0.414 day orbital 
period; however this mainly refiects the contribution of the third component because spectral lines 
of the contact binary are highly broadened by the fast rotation and are practically invisible in 
classification dispersion spectra. The 2MASS infrared color, J — K = 0.282, corresponds to the 
F5V spectral type and our own classification spectra also indicate the F5V type. The Hipparcos 
parallax of the system, 5.16 =b 1.69 mas is too small and uncertain to draw any conclusions. 



2.4. VZ CVn 



VZ CVn (HD 117777, HIP 6601 7) is a rather bright 



= 9.35), detached eclipsing binary. 
It was discovered to be variable by IStrohmeier fc Kniggd (jl960l ). The system is k r iown to show 
a slow intrinsic variability in the LC (see |Poppeiill988l ). Recently, llbanoglu et al.l (j2007l ) - after 
removing the eclipse variability - detected 7 Dor-type oscillations of the primary component and 
determined the first set of photometric elements. As this star has an orbital period of 0.842 days, 
it is difficult to cover all orbital phases for this target from one site; in addition intrinsic night-to- 
night variations of a few hundreds of a magnitude make a consistent LC solution almost impossible 
using ground-based photometry. The system is therefore a prime candidate for satellite, continuous 
photometric observations. 



VZ CVn was observed spectroscopically by iPoppeij (|l988l ). who determined the first radial 
velocity orbit from photographic spectra: Ki = 144.1 ± 1.2 km s~^, K2 = 185.4 it 3.0 km s~^ 
and the separate systemic velocities of the two components as Vqi = —20.7 it 0.7 km s~^, and 
V02 = -21.7 ±2.0 km s'^ 

Our observations provide v ery well define d BFs leading to spectroscopic elements that are 
mostly consistent with those of iPoppeii (jl988l ) . The largest discrepancy is the smaller value of 
the semi-amplitude K2. The spectroscopic orbit is, however, of excellent quality, leading to a 
total projected mass of (Mi + M2)sin^i = 2.676 it O.OO7M0, which corresponds to the relative 
precision as good as 0.26 %. The broadening functions are well separated in orbital quadratures, 
thus supporting the photometric classification of the system as a detached binary. 

The 2MASS infrared color of the system, J — K = 0.216 indicates the F2V spectral type. Our 
classification spectra correspond to FOV. The Hipparcos parallax is small, 2.23 it 1.22 mas and too 
inaccurate to be of use. 
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2.5. GK Cep 



The bright variable star GK Cep (HD 205372, HIP 106226; Vmax = 6.99), originally classified 
as an RR Lyr p ul sator, was later found to be an eclipsing binary with a 0.483 days period (see 
Strohmeieiill963l ). iBartollini et al.l (|l965l ) observed the system photometrically and spectroscopi- 
cally and determined the correct orbital period, P = 0.936171 days, and the first set of spectroscopic 
elements: Vq = —22 km s~^, Ki = 172 km s~^, K2 = 187 km s~^. The authors classified GK Cep 
as a /3 Lyrae variable and, with the photometrically determined inclination angle, i = 71°, they de- 
termin ed the masses of the components, Mi = 2.7 Mq and Mi = 2.5 Mq. Later. iHutchings &: Hill 



(jl973l ) analyzed the LCs of the system assuming the Roche model and found that the system is a 
close, but detached binary; the more massive component is the cooler one. 



The orbital period analysis of iKreiner et al.l (|l99d ) suggests the presence of an invisible third 
body in the system on a 18.8 year orbit. Assuming that the orbits are co-planar, the authors 
estimated the mass of the third component as 1.34 Mq. 

Our spectroscopy confirms the existence of the third component. It is, however, faint with its 
light contribution as small as L^/{Li + L2) ~ 0.025. The profile of the third component is not well 
visible in all of the BFs because it blends with the more massive component, especially around the 
second quadrature. This circumstance complicated our usual approach, i.e. the fitting of a triple 
Gaussian model to the observed BFs. Hence we averaged all BFs (in the heliocentric RV system) 
and the mean BF was used to define the third component contribution which was subsequently 
subtracted from all individual BFs. The RV of the third component, RV^ = —19.3 km s~^, is rather 
different from the systemic velocity of the system, Vq = —31.06 it 0.45 km s~^, indicating that the 
faint companion is not physically related to the binary. 



Our spectroscopic elements (Table [2]), largely agree with those of IBartollini et al.l (|l965l ). Our 
RV semi-amplitudes are, however, smaller re sulting in an est ir nate o f the total mass of the system 
that is smaller by about 15% than that of IBartollini et al.l (119651 ). The very large mass ratio, 
q = 0.913 =fc 0.005, wh ich would be unusual if it were a contact binary, supports the result of 
Hutchings Hilll (119731 ) that the system is detached. 



The Hipparcos parallax vr = 5.17 it 0.32 mas (the distance d = 193 it 11 pc) is precise thanks 
to the high ecliptical latitude. The 2MASS infrared color, J — K = 0.007 corresponds to the AlV 
spectral type, which is in good accord with the previous spectral classification, A2V (Hill et al., 
1975). Our classification spectra give a slightly earlier spectral type, AOV. If the components of 
GK Cep are normal main-sequence stars, the best accord with the Hipparcos parallax is achieved 
for the A2V spectral type. 
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2.6. RW Com 



RW Com (HIP 61243, Vmax ~ H-O mag, sp. type G5 - G8) is one of the s h ortest -period 
(P = 0.237346 day) W UMa systems. Its variabihty was first noticed by iJordanI (jl923l ). This 
late- type contact binar y is known to show an enhanced surface activity resulting in an asymmetric 



LC ( Milone et al.lll987 



The orbi tal period of the system shows a continuous decrease of dP/dt 



6.06 X 10 * day/year ( Qianll2002l ) accompanied by a sinusoidal variation with a 16-year periodicity. 



The first spectroscopy of R W Com was carried out by IStruvd (jl95d ) who found the H & K 
Ca II lines in emission. Later, iMilone et al.l (|l985l ) presented the first RV orbit of the system: 



K1 + K2 = 304±5 km s'^ , {Mi + M2)sui^i = O.69±O.O4M0, M1/M2 = 2.9±0.2 (this corresponds 
io q = 0.349 ± 0.022) and Vq = -53 ± 4 km 8"^ Unfortunately, the spectra were of a poor quality, 
with relatively long exposure times typically lasting 20 minutes, which resulted in heavy blending 
of the cross-correlation functions used for radial velocity determinations. 

Because of the short orbital period of the system, we kept all our exposures eq ual to 500 sec 
(2.43% of the orbital period). Our data, however, by far supersede those presented by 



1985 ). The resulting parameters given in Table [2] are inconsistent with those given by 



Milone et al 



Milone et al 



19851 ): the most striking difference is our much higher mass ratio, q = 0.471 it 0.006, and the 



substantially larger total projected mass, (Ml + M2) sin^ i = 1.052 ± O.O13M0. It is interesting to 
note, that very similar mass ratios, all close to 0.5, have been found for tw o other contact, very 
late-type binaries with extrem ely short orbital periods: C C Com {q = 0.527, iPribulla et al.l 120071 ) 
and GSC 1387 475 {q = 0.474, iRucinski Pribullalbood ). The broadening functions do not show 
any trace of the third component which was indicated by the cyclic period change. This sets the 
upper limit on its luminosity at about L^/{Li + L2) < 0.03. 

Using the Hipparcos parallax, tt = 14.33 it 3.36 mas, and the know n maximum apparent 
magnitude, Vmax = 11-0, one obtains My = 6.8±0.8. With the calibration of lRucinski &: Duerbeck 



(119971) and using 5-1/ = 1.07 from the TYCH02 catalogue, one obtains M^""^ = 6.12. The 2MASS 



infrared color of RW Com, J—K = 0.618 corr esponds to the K2V spectral type, which is inconsistent 
with the G5-G8 spectral type estimated by iMilone et al.l (|l985l ). It is however, in better accord 
with the very short orbital period of the binary and our rather rough estimate of K2/5V. 



2.7. V2610 Oph 



V2610 Oph (HD 162905, SAO 141948) was found to be variable bv lWils k DworakI (|2003l ) using 
the Stardial images. The authors classified it as a W UMa-type binary with an amplitude of about 
0.16 mag and noted that its spectral type, KO, is too late for the observed orbital period. Later 



Tas &: EvrenI (120061 ) obtained precise BVR light curves and determined the preliminary geometric 
parameters of the system as q = 0.55, i = 54.22° and found that the system is detached but close 
to contact. The authors assumed no third light in the system. No spectroscopic observations of 
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the system have been pubhshed yet. 



Our spectroscopy immediately revealed a rather different pictu re: V2610 Oph is a quadruple 
system. It is somewhat similar to VW LMi (jPribulla et al.l 120061 ) as it consists of an eclipsing 
contact binary and a detached non-eclipsing pair. The orbit of the second binary, with a period of 
P = 8.47 days, is slightly eccentric (see Table H] and Fig. [5]). 



The multiplicity of V2610 Oph makes the photometric solution of iTas EvrenI (|2006l ) com- 
pletely inapplicable. The spectroscopic mass ratio, q = 0.289 and the large projected masses of the 
components, (Mi + M2) sin^ i = 1.408Mq, indicate a much higher inclination angle than what was 
found previously. The light contribution of the second binary in V2610 Oph around the maxima 
of the eclipsing pair, as found by modeling the BFs, is (L3 + Li)/{Li + L2) = 1.28 =b 0.11, i.e., the 
non-eclipsing pair is brighter than the contact b inary. When we correct the observed photometric 
amplitude, Al/ = 0.163 mag (|Tas &: Evrenll2006l ) for this contribution, we obtain the full amplitude 
of the binary as 0.41 it 0.02 mag, indicating a high orbital inclination angle. 

The center-of-mass RV of the second pair, Vq = 60.76 it 0.11 km s^^ is fairly close to the 
systemic velocity of the eclipsing binary, 71.70 it 1.27 km s~^. The RV difference indicates mutual 
orbital motion of the binaries. No systemic- velocity changes in either of the binaries were noted 
during our observing run that extended for nearly one year. Hence, the outer orbital period is at 
least several years long. The angular separation between the binaries must be smaller than about 
1 arcsec, because V2610 Oph appeared to be a single star on our spectrograph slit even during 



the n ights of excellent seeing; besides, the system is not known to be a visual binary (iMason et al, 
200lh. 



A new LC analysis using the new spectroscopic mass ratio and including the third light is 
necessary to determine the inclination angle. Then, having the masses of the close binary, masses 
of all four compone nts could be determined using the "linked mass-ratio" technique as used in 
Pribulla et al] (|20od l for VW LMi. 



The infrared color of the system, J — K = 0.393, corresponds to G4V spectral type, while the 
TYCHO color, B — V = 0.587 corresponds to GOV spectral type. A direct spectral classification is 
difficult to perform but gives a reasonable range of admissible types is F8 - G2V which is consistent 
with the projected mass of the dominant primary component of the contact binary. Mi = LOOM©. 
V2610 Oph was not observed by the Hipparcos satellite and its trigonometric parallax is unknown. 



2.8. V1387 Ori 

V1387 Ori (HD 42969, HIP 29186) is a member of a trapezoidal, visual, multiple system. 
The other members are GSC 1318 59 and GSC 1318 317; the latter star forms the close visual 
binary HDS 838. Neither of these visual companions entered our spectrograph slit. Establishing 
the physical bond between the members of this mini-cluster would require precise proper motions 
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and RVs. 



The variability of V1387 Ori was discovered during the Hipparc os mission (ESA 



19971). In the 



Hipparcos catalogue it was classified as a /3 Lyrae system, but later iDuerbeckl (jl997l ) suggested it 
to be a contact, W UMa-type system. The Hipparcos light curve shows a maximum following the 
primary minimum that is substantially brighter than the other maximum (similar to the contact 
binaries AG Vir and DU Boo). No ground-based photometric nor spectroscopic observations of the 
system have been published yet. 

Our spectroscopy revealed that V1387 Ori is a triple system. The radial velocity of the third 
component appears to vary as indicated by our first spectra obtained at HJD 2454073.7873, giving 

= 34.43 km s~^. This observation was well separated in time from the remaining observations 
which give a constant velocity of -RV3 = 24.4 ± 0.4 km s^^. The contribution of third component is 
affected by the large O'Connell effect observed in the system: at Max I, L3/(Li +L2) = 0.128±0.012 
while at Max II, L^/{Li + L2) = 0.186 it 0.023. When the third-body signatures are removed from 
the BFs, the orbit of the eclipsing pair is well defined with a low mass ratio, q = 0.165 it 0.005, 
which would be rather typical for an A-type contact binary. No surface inhomogeneities on either 
of components were noticed in the BFs in spite of the large O'Connell effect. 

The Hipparcos parallax for this system is negative, vr = — 7.86±1.96 mas, which most probably 
results from V1387 Ori being located in the visual multiple system. The 2MASS color, J — K = 
0.155, corresponding to A2V is consistent with our spectral classification, A4V. 



2.9. AU Ser 



AU Se r (GSC 1502 1762) is a rather faint iVmax = 10.9) con tact binary discov ered by 
Hoffmeisterl (j 19351 ). The first photoelectric photometry presented by iBinnendijkl (j 19721 ) showed 
a significant asymmetry in the LC (Max I - Max II = —0.05 mag ) and a 0.2 mag difference in 
the dept h of th e eclipses. The light curves obtained by IBinnendijkl (jl972l ) were later analyzed by 
Kaluznyl (jl986l ) who interpreted the large asymmetry as a hot spot located close to the "neck" 
connecting the components. The best solution, under the assumption of the A-type of the light 
curve, was obtained for a rather large mass ratio q = 0.80. The authors, however, admitted tha t 
the acceptable range for the mass ratio is rather wide, between 0.70 and 1.15. Later iGuroll (|2005l ) 
analyzed photometric observations of this target covering the period from 1969 until 2003. Their 
analysis of the observed times of minima indicates the presence of a third body on a 94.15 year orbit 
with an estimated minimum mass of = 0.53Mq. According to the authors the differences in 
the LC maxima levels appear to be cyclic on a time scale of about 30 years. 



Spectroscopic observations of this system were performed by iHrivnakI (|l993l ). who reported 



q = 0.71 and a projected total mass of the system of (Mi + M2)siv?i = I.SIMq. The systemic 
velocity of the system was not given. Our new spectros copic obs e rvatio n leads to a projected mass 
and a mass ratio within the errors of those found by iHrivnakI (jl993l ). We do not see evidence 
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for a third component in the BFs. The phased BFs, however, show a large cool spot on the 
secondary component which is visible after the second quadrature (see Fig. H|). This dark, localized 
spot significantly deforms the BFs and affects the determined RVs. A reliable determination of 
the spectroscopic elements would require modeling of simultaneous spectroscopic and photometric 
observations. 

The system appears to be of relatively late spectral type, G4V, for its orbital period of 0.386 
d. AU Ser was too faint for the Hipparcos mission, hence its parallax is unknown. 



2.10. FT UMa 



FT UMa (HD 75840, HIP 43738) was discovered to be v ariable by the Hipparco s satellite. 
In the special General Catalogue of Variable Stars namelist (IKazarovets et al.lll999l ). the RRc 
classification was suggested with a pulsational period of 0.3273519 days . High precision BVR 
photometric observations of FT UMa were published by ISelam et al.l (|2007l ) . The authors correctly 
classified the system as a close binary, obtained a low mass ratio of q = 0.25 ib 0.01, a low orbital 
inclination of z = 54.48 ± 0.80°, and suggested a contact configuration for the system. In spite of 
the fact that this target is relatively bright, Vmax = 9.25, and the photometric amplitude is as large 
as 0.17 mag, no other observations of FT UMa have been published yet. 

Our spectroscopy clearly shows that FT UMa i s a triple system containing a close eclipsing 
binary with a period of P = 0.6547038 days. FT UMa is not listed in the WDS catalogue as a visual 
binary. During our spectroscopic observations the star always appeared as a single object, which 
sets an upper limit to the angular separation of the components at about 1 arcsec. The center- 
of-mass velocity of the close binary, Vq = —33.66 it 1.30 km s~^, is significantly different from the 
third-component velocity which slowly changed from —22 km s~^ to about —10 km s~^ during our 
run. Unfortunately, the RVs of the close pair are rather imprecise to resolve the hierarchy and 
multiplicity of the system: it can either be a hierarchical triple or a quadruple system consisting 
of two binaries (with the non-eclipsing pair being the SBl). The light contribution of the third 
component around the maxima of the eclipsing pair is L^/[Li + L2) = 1.01 ± 0.03. The full, 
undiluted, photometric amplitude of the eclipsing binary would be then about 0.37 mag. The 
presence of such a strong third light made the photometric solution of (ISelam et al.ll2007l ) of limited 



use. 

The profiles of all three components are well separated in the extracted BFs, indicating that 
the close eclipsing binary is very probably detached, which is consistent with the large mass ratio, 
q = M2/M1 = 0.984 it 0.019, but rather atypical for contact binaries. On the other hand, Hipparcos 
photometry phased with the adopted period of 0.6547038 days results in a LC that is typical of 
contact binaries. 



The TYCH02 color of the system, B — V = 0.404 reflects the combined contributions of the 
eclipsing pair and of the third component. Our spectral classification spectra suggest the FOV 
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spectral type, which is consistent with the 2MASS J — K = 0.207. The Hipparcos parallax is fairly 
small, 6.90 it 1.28 mas, and too inaccurate to be of use. 



3. SUMMARY 

With the ten new short-period binaries, this paper brings the number of the systems studied 
at the David Dunlap Observatory to one-hundred thirty. With the closure of the observatory, this 
series is coming to an end. However, the number of known, bright close binaries that remained 
unobserved at DDO is moderate; only a dozen or so known W UMa-type eclipsing binaries brighter 
than about V = that were accessible from DDO remained. We plan to publish the data for the 
unfinished cases in the last, fifteenth installment of this series. 

The highlights of the current series are: (1) the discoveries of two quadruple systems TZ Boo 
and V2610 Oph, (2) the discoveries of four triple systems: EL Boo, GK Cep, V1387 Ori, and 
FT UMa, (3) the spotted contact binary AU Ser, with a large spot on the secondary component. 
Non e of the dete c ted sp ectroscopic multiple systems were previously noted to be visual binaries 



(see iMasonetaLl (1200l|)). 



While for four systems, EL Boo, V2610 Oph, V1387 Ori and FT UMa, we are presenting 
the first spectroscopic observations, the quality of the data for the remaining six systems is much 
improved relative to the previously published investigations. 

Numerous discoveries and reliable solutions of triple and quadruple systems show that the BE 
deconvolution approach utilizing the SVD method is a powerful technique. In this series, good 
examples are the three triple systems with a dominant third component, EL Boo, V2610 Oph and 
FT UMa. Also, the BE technique enabled us to determine the first reliable set of spectroscopic 
elements for TZ Boo and to reveal the true nature of this system. The third component of this 
system, moving in 9.48 day orbit, blends with the components of the close binary, a condition 
that made all previous spectroscopic studies of TZ Boo unsuccessful. The second binary may be 
the cause of the large perturbations observed in the light curve of the eclipsing pair. However, 
with several multiple systems analyzed in the DDO papers, we also see the weakness of t he B E 



technique that should be addressed in the future: As described in iRucinski &: Pribullal (|2008l ). 
the contribution of the third component to the total light of the system is always overestimated. 
There are two reasons for this discrepancy: (i) the continuum spectrum rectification differently 
affects spectra of (usually) sharp-lined companions and of heavily broadened short-period binaries; 
(ii) the technique must use a single template spectrum, but this affects the luminosity ratio when 
the spectral types of the third component and the contact binary are significantly different. 
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Captions to figures: 

Fig. 1. — Radial velocities of the systems TZ Boo, VW Boo, EL Boo, and VZ CVn are plotted in 
individual panels versus the orbital phases. The lines give the respective circular-orbit (sine-curve) 
fits to the RVs. TZ Boo is a quadruple system consisting of a contact and a detached single-line 
binary while VW Boo is a contact binary. EL Boo is a triple system harboring a close binary and 
VZ CVn is a detached or a semi-detached binary. The circles and triangles in this and the next 
two figures correspond to components with velocities Vi and V2, as listed in Tabled! respectively. 
The component eclipsed at the minimum corresponding to Tg (as given in Table [5]) is the one that 
shows negative velocities for the phase interval 0.0 — 0.5 and is the more massive one. Short marks 
in the lower parts of the panels show phases of available observations which were not used in the 
solutions because of the spectral line blending. 

Fig. 2.— The same as for Figure E but for GK Cep, RW Com, V2610 Oph, and V1387 Ori. 
GK Cep is a close binary with a faint third component; RW Com is a contact binary; V2610 Oph is 
a quadruple system consisting of a contact eclipsing binary and a wide non-eclipsing pair; V1387 Ori 
is a triple system containing a close eclipsing binary. V1387 Ori forms a part of a trapezoidal visual 
multiple system. 

Fig. 3. — The same as for Figures [1] and [21 for the two remaining systems AU Ser, and FT UMa. 
AU Ser is a contact binary while FT UMa is a close eclipsing binary in a triple system. 

Fig. 4. — The broadening functions (BFs) for all ten systems of this group, selected for orbital phases 
close to 0.25 or 0.75. The phases are marked by numbers in the individual panels. Additional 
components to the close binaries, TZ Boo, EL Boo, V2610 Oph, V1387 Ori, and FT UMa, are 
strong and clearly visible. The third component in GK Cep is a small peak close to the binary 
center-of-mass velocity. All panels have the same horizontal range, —500 to +500 km s~^. 

Fig. 5. — The radial velocities and the best fits for the third component of TZ Boo and the second 
non-eclipsing binary in the quadruple system V2610 Oph. 
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Fig. 1. — Radial velocities of the systems TZ Boo, VW Boo, EL Boo, and VZ CVn are plotted in 
individual panels versus the orbital phases. The lines give the respective circular-orbit (sine-curve) 
fits to the RVs. TZ Boo is a quadruple system consisting of a contact and a detached single-line 
binary while VW Boo is a contact binary. EL Boo is a triple system harboring a close binary and 
VZ CVn is a detached or a semi-detached binary. The circles and triangles in this and the next 
two figures correspond to components with velocities Vi and V2, as listed in Tabled! respectively. 
The component eclipsed at the minimum corresponding to Tq (as given in Table [2]) is the one that 
shows negative velocities for the phase interval 0.0 — 0.5 and is the more massive one. Short marks 
in the lower parts of the panels show phases of available observations which were not used in the 
solutions because of the spectral line blending. 
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Fig. 2.— The same as for Figure [H but for GK Cep, RW Com, V2610 Oph, and V1387 Ori. 
GK Cep is a close binary with a faint third component; RW Com is a contact binary; V2610 Oph is 
a quadruple system consisting of a contact eclipsing binary and a wide non-eclipsing pair; V1387 Ori 
is a triple system containing a close eclipsing binary. V1387 Ori forms a part of a trapezoidal visual 
multiple system. 
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Fig. 3. — The same as for Figures [T] and [21 for the two remaining systems AU Ser, and FT UMa. 
AU Ser is a contact binary while FT UMa is a close eclipsing binary in a triple system. 
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Fig. 4. — The broadening functions (BFs) for all ten systems of this group, selected for orbital phases 
close to 0.25 or 0.75. The phases are marked by numbers in the individual panels. Additional 
components to the close binaries, TZ Boo, EL Boo, V2610 Oph, V1387 Ori, and FT UMa, are 
strong and clearly visible. The third component in GK Ccp is a small peak close to the binary 
center-of-mass velocity. All panels have the same horizontal range, —500 to +500 km s~^. 
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Fig. 5. — The radial velocities and the best fits for the third component of TZ Boo and the second 
non-eclipsing binary in the quadruple system V2610 Oph. 
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Table 1. DDO radial velocity observations (the full table is available only in the electronic form) 



Tari 


get 


HJD-2,400,000 


Vi 


Wi 


V2 


W2 


Phase 








[km s"^] 




[km s-i] 






VW 


Boo 


54558.8318 


0.00 


0.00 


0.00 


0.00 


0.9153 


VW 


Boo 


54558.8472 


0.00 


0.00 


0.00 


0.00 


0.9602 


VW 


Boo 


54558.9127 


-68.22 


0.50 


216.48 


0.50 


0.1516 


VW 


Boo 


54558.9233 


-67.24 


1.00 


237.55 


1.00 


0.1826 


VW 


Boo 


54559.7940 


116.90 


1.00 


-200.04 


1.00 


0.7261 


VW 


Boo 


54559.8047 


116.91 


1.00 


-208.79 


1.00 


0.7574 


VW 


Boo 


54563.6869 


-51.23 


1.00 


165.41 


1.00 


0.0984 


VW 


Boo 


54563.6977 


-57.31 


1.00 


197.76 


1.00 


0.1299 


VW 


Boo 


54563.7091 


-64.92 


1.00 


235.53 


1.00 


0.1632 


VW 


Boo 


54563.7200 


-71.21 


1.00 


243.67 


1.00 


0.1951 



Note. — The table gives the RVs Vi for observations described in the paper. 
The first 10 rows of the tabic arc shown for one of our typical program stars, 
VW Boo. For first program star, TZ Boo, where phase smoothed BFs were used, 
the heliocentric Julian dates are not given. Observations leading to entirely insep- 
arable broadening function peaks arc given zero weight; these observations may 
be eventually used in more extensive modeling of the broadening functions. Zero 
weights were assigned to observations of marginally visible peaks of the secondary 
(sometimes even the primary) component. The RVs designated as V\ correspond 
to the more mgissive component; it was always the component eclipsed during the 
minimum at the epoch To (this does not always correspond to the deeper mini- 
mum and photometric phase 0.0). The phases correspond to spectroscopic To and 
periods given in Table 2, but not necessarily to the photometric ephemerides given 
below the table). 
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Table 2. Spectroscopic orbital elements 



Name 


Type 


Other names 


Vb 


Ki 


£1 


To - 2,400,000 


P (days) 


9 




Sp. type 






K2 


62 


(C - G)(d) [tj 


(Ml + M'2 ) sm I 




TZ Boo 


EW(A/W) 


BD+40 2857 


-46.57(0.90) 


57.8(1.45) 


8.70 


54042.7482(4) 


0.2971597 


0.207(5) 




F/G5 


HIP 74061 




280.02(1.50) 


10.49 


+0.0002 [+5.191.0] 


1.188(18) 




VW Boo 


EW 




22.79(0.66) 


101.39(1.07) 


6.65 


54573.9227(4) 


0.3423157 


0.428(5) 




G5V 


HIP 69826 




236.74(1.0.3) 


5.99 


-0.0032 [+6,058.5] 


1.371(14) 




EL Boo 


EW 


BD+14 2788 


-24.58(1.12) 


64.19(1.70) 


9.58 


54584.1779(9) 


0.413772 


0.248(7) 




F5V 


HIP 72391 




259.07(1.89) 


9.31 


-0.0850 [+14,704] 


1.448(28) 




VZ CVn 


EB 


HD 117777 


-22.99(0.16) 


141.51(0.26) 


1.64 


54573.9540(3) 


0.84246123 


0.8252(19) 




FOV 


HIP 66017 




171.48(0.25) 


1.72 


+0.0009 [+2,461] 


2.676(7) 




GK Cep 


EB 


HD 205372 


-31.06(0.45) 


163.44(0.73) 


5.05 


54464.4014(8) 


0.936169 


0.913(5) 




AOV 


HIP 106226 




178.95(0.73) 


4.20 


-0.0062 [+2,098] 


3.893(26) 




RW Com 


EW(W) 




39.35(0.86) 


112.04(1.28) 


6.33 


54272.7645(3) 


0.2373464(7) 


0.471(6) 




K2/5V 


HIP 61243 




237.70(1.29) 


7.39 


+0.0005 [+7,465.5] 


1.052(13) 




V2610 Oph 


EW{A) 


HD 162905 


71.70(1.27) 


72.06(2.06) 


11.87 


54461.5659(10) 


0.426512(3) 


0.291(9) 




F8/G2V 






247.42(2.08) 


12.11 


+0.0109 [+4,904] 


1.441(31) 




V1387 Ori 


EW? 


HD 42069 


56.92(0.65) 


33.21(0.98) 


4.37 


54204.0187(13) 


0.730166 


0.165(5) 




A4V 


HIP 29186 




200.72(1.18) 


8.50 


+0.0082 [+7,811] 


0.969(16) 




AU Ser 


EW(A) 




-61.59(0.84) 


138.77(1.33) 


6.36 


54597.4761(7) 


0.386499 


0.709(8) 




G4V 






195.64(1.34) 


10.43 


+0.0042 [+5,426] 


1.498(19) 




FT UMa 


EB 


HD 75840 


-33.66(1.30) 


155.15(2.14) 


10.78 


54561.9614(19) 


0.654704 


0.984(19) 




FOV 


HIP 43738 




157.73(2.13) 


10.88 


-0.0136 [+9,258.5] 


2.077(45) 





Note. — The spectral types given in column 2 relate to the combined spectral type of all components in a system; they are given in 
parentheses if taken from the literature, otherwise they are new. The convention of naming the binary components in the table is that the 
more massive star is marked by the subscript "1", so that the mass ratio is defined to be always q < 1- The standard errors of the circular 
solutions in the table are expressed in units of the last decimal places quoted; they are given in parentheses after each value. The center-of- 
mass velocities (Vq), the velocity amplitudes {Ki) and the standard unit-weight errors of the solutions (e) arc all expressed in km s^-'^. The 
spectroscopically determined moments of primary or secondary minima are given by Th (corresponding approximately to the average Julian 
date of the run); the corresponding (O — C) deviations (in days) have been calculated from the available prediction on primary minimum, 
as given in the text, using the assumed periods and the number of epochs given by [E]. The values of (Mi + M2) sin^ i are in the solar mass 
units. 

Ephemerides {HJDmin ~ 2,400,000 + period in days) used for the computation of the (O — C) residuals: 
TZ Boo: 52500.1920 + 0.2971597 

VW Boo: 52500.0062 + 0.3423157 
EL Boo: 48500.1594 + 0.413772 
VZ CVn: 52500.6560 + 0.84246123 
GK Cep: 52500.325 + 0.936169 
RW Com: 52500.1432 + 0.2373463 
V2610 Oph: 52369.95 + 0.42651 
V1387 Ori: 48500.6839 + 0.730166 
AU Ser: 52500.3392 + 0.386497 
FT UMa: 48500.3999 + 0.6547038 
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Table 3. Radial velocity observations (the full table is available only in electronic form) of third 

and fourth components of multiple systems 



Target 


HJD-2,400,000 


Vs 


V4 






[km s-i] 


[km s"^] 


V2610 Oph 


54302.7053 


91.91 


27.54 


V2610 Oph 


54302.7160 


91.99 


28.35 


V2610 Oph 


54306.6868 


17.32 


107.57 


V2610 Oph 


54306.6976 


16.85 


106.22 


V2610 Oph 


54306.7090 


20.25 


108.84 


V2610 Oph 


54306.7197 


17.23 


106.02 


V2610 Oph 


54307.5974 


48.53 


76.07 


V2610 Oph 


54307.6080 


48.49 


75.25 


V2610 Oph 


54307.6193 


49.66 


75.37 


V2610 Oph 


54307.6299 


49.44 


73.04 



Note. — The table gives the RVs Vi for the third and 
fourth components. The typical 10 rows of the table for 
quadruple system, V2610 Oph, are shown. Observations of 
the quadruple system V2610 Oph leading to entirely insep- 
arable broadening function peaks of the components of the 
second binary have been omitted from the table and are 
not used in the computation of the orbit. 



Table 4. Spectroscopic orbital elements of the second non-eclipsing binaries in the quadruple 
systems TZ Boo and V2610 Oph . Orbit of the second binary in TZ Boo is circular, thus 

634 = 0.00 and 0^34 = 7r/2 



Parameter 



TZ Boo 



P34 [days] 


9.47765 


0.00034 


To [HJD] 


2 454 042.818 


0.007 


Vo [km s-i] 


-54.64 


0.12 


K:i [km s-i] 


43.15 


0.14 


03 sini [Rq] 


8.08 


0.03 


/(m) [Mo] 


0.0793 


0.0008 


/2610 Oph 






P34 [days] 


8.47093 


0.00025 


e34 


0.073 


0.003 


CJ [rad] 


5.88 


0.04 


To [HJD] 


2 454 461.85 


0.06 


Vo [km s"'] 


60.76 


0.11 


K3 [km g-^j 


57.77 


0.23 


K4 [km s"^] 


62.15 


0.24 


g = ^^4/7^3 


0.929 


0.005 


(03 + 04) sini [R0] 


20.01 


0.06 


(Ms + M4) sin^ i [Mq] 


1.502 


0.013 



Note. — The table gives spectroscopic elements 

of the second binaries in TZ Boo and V2610 Oph: 
orbital period (-P34), eccentricity (634), longitude 
of the periastron passage (w), time of the peri- 
astron passage (To), systemic velocity (Vo), semi- 
amplitudes of the RV changes (if 3,1^4). Cor- 
responding mass ratio q, and projected relative 
semi-major ((03 + 04) sini) and total mass ((M3 -|- 
M4) sin^ i) is given for V2610 Oph where both lines 
of the second binary could be measured. For single- 
lined non-eclipsing binary in TZ Boo only 03 sin i 
and f{m) is given. 



